Entanglement creation in circuit QED via Landau-Zener sweeps 
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A qubit may undergo Landau-Zener transitions due to its coupling to one or several quantum 
harmonic oscillators. We show that for a qubit coupled to one oscillator, Landau-Zener transitions 
can be used for single-photon generation and for the controllable creation of qubit-oscillator en- 
tanglement, with state-of-the-art circuit QED as a promising realization. Moreover, for a qubit 
coupled to two cavities, we show that Landau-Zener sweeps of the qubit are well suited for the 
robust creation of entangled cavity states, in particular symmetric Bell states, with the qubit act- 
ing as the entanglement mediator. At the heart of our proposals lies the calculation of the exact 
Landau-Zener transition probability for the qubit, by summing all orders of the corresponding series 
in time-dependent perturbation theory. This transition probability emerges to be independent of 
the oscillator frequencies, both inside and outside the regime where a rotating-wave approximation 
is valid. 
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I. INTRODUCTION 

Entanglement is a purely quantum mechanical prop- 
erty of multipartite systems. A system is entangled if 
its quantum state cannot be described as a direct prod- 
uct of states of its subsystems. Entanglement is mea- 
surable in terms of nonclassical correlations of the sub- 
systems. Famous examples are the Einstein-Podolsky- 
Rosen correlations between positions and momenta of 
two particles and the violation of the Bell inequal- 
ities by spin systems that are described by Bell states 
• Fundamental tests of entanglement and nonlocality 
were performed in quantum optics [3|, with the intrigu- 
ing outcome that measured nonclassical correlations be- 
tween entangled spatially separated subsystems rule out 
local realism [4|. 

With the advent of quantum information theory in re- 
cent years, the interest in entanglement has broadened. 
Many efforts exist to make use of entanglement in in- 
formation processing or to quantum-communicate with 
built-in security With quantum information pro- 

cessing in mind, in this paper we propose to create en- 
tanglement in two spatially separated elements (circuit 
cavities) of a superconducting circuit by letting a third 
element (the superconducting qubit) undergo a Landau- 
Zener sweep. This will be a robust method to create 
Bell states in two-cavity circuit cavity quantum electro- 
dynamics (QED). 

In optical cavity QED, atoms (qubits or N- level sys- 
tems) become entangled with optical cavity modes (quan- 
tum harmonic oscillators). The creation of atom-cavity 
entanglement in cavity QED is possible because the 
strong coupling regime can be realized, where Rabi oscil- 
lations between atomic and optical excitations occur on 
a faster time scale than spontaneous emission and cav- 
ity decay. By adiabatic passage the quantum state of 



an atom flying through an optical cavity can be mapped 
onto the quantum state of the cavity. It was investigated 
theoretically @, 0] and shown experimentally Q that in 
a similar manner also two modes of the same cavity can 
be entangled. Alternatively, two spatially separated op- 
tical cavities could be entangled by letting an atom fly 
successively through both cavities [§, [l(| , or by detecting 
photons leaving the two cavities . 

If the aim is to build devices for quantum information 
processing, then scalability to several qubits and oscil- 
lators is an important prerequisite. It seems technolog- 
ically very challenging to scale up optical cavity QED. 
Recently, a new research field called circuit QED has 
emerged in which analogues of cavity QED have been 
realized with superconducting qubits [H, [lU and os- 
cillators. A flux qubit was coupled strongly to a su- 
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perconducting quantum interference device 
charge qubit to a transmission line resonator 
and Rabi oscillations have been observed in both exper- 
iments. There are analogous proposals to couple super- 
conducting qubits to nanomechanical resonators (l7l Il8| . 
Superconducting circuits are promising because of their 
potential scalability and because many of their parame- 
ters are tunable over a broad range. 

Scaling up present-day circuit QED can be done by 
increasing the number of qubits, the number of oscilla- 
tors, or both. The first option is encountered in most 
proposals, where the quantum harmonic oscillator is to 
be used as a bus to manipulate and read out the qubits 
[HI GJ, GJ, S3 • The second option is the circuit analogue 
of two-cavity (or AT-cavity) QED [i| |T(| , and the present 
paper belongs to this category. Also in this group is the 
recent proposal to couple a superconducting qubit to two 
transmission-line resonators with orthogonal field polar- 
izations plj . As to the third option, examples of propos- 
als for experiments with several superconductin g q ubits 
coupled to several oscillators can be found in [LSI . S3- 
Probably the first realization of such a complex setup 
will consist of two qubits independently coupled to their 
own oscillators, allowing independent qubit state manip- 
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ulation and readout. 

One method to manipulate the state of an isolated 
qubit is to use Landau- Zener sweeps [H, [3, HH, as 
discussed in more detail below. LZ transitions can be 
used to control qubit gate operations [H, Ht| and to 
read out qubits [281 ] - Recently, LZ transitions have been 
observed in various experiments with superconducting 

qubits minim m. 

In this paper we concentrate on quantum state ma- 
nipulation in multi-cavity circuit QED via Landau-Zcncr 
sweeps of a qubit. Bit flips in the qubit can take place 
even in the absence of a direct coupling between the qubit 
levels, induced instead by the coupling to the oscilla- 
tors. Oscillators in highly excited coherent states can 
be described classically and give rise to a non-monotonic 
LZ transition probability as a function of the coupling 
strength [33[. Here we focus on the experimentally rel- 
evant situation in which the oscillators all start in their 
ground states. We show how single photons can be cre- 
ated, not only in the presence of only one oscillator as in 
our prior work [34], and how symmetric Bell states can 
be created in two circuit oscillators. 

The decisive advantage of our proposal is that qubit- 
oscillator interaction strengths are static, in contrast to 
the situation in standard cavity QED where these inter- 
actions require precise dynamical control [llSHcl- Other 
advantages of LZ transitions are that even adiabatic in- 
teractions can be performed rather fast 35|, and that 
LZ transition probabilities are extremely robust under 
dephasing |3a.l37|. 

In Sec. [IT] we review the standard LZ problem, before 
deriving in Sec. IIIII exact LZ transition probabilities for 
a qubit coupled to N harmonic oscillators. This result 
is used in Sec. IIVI about single-photon generation and in 
the central section [V] about the generation of entangled 
cavity states via LZ sweeps of the qubit. The realization 
in circuit QED of our proposals is discussed in Sec. IVI1 
before we conclude. 



II. STANDARD LANDAU-ZENER TRANSITION 

For later use, we shall briefly review the well-known 
Landau-Zener transition in a qubit. Consider a two-state 
system with states |f) and |J.). The standard Landau- 
Zener (- Stueckelberg) problem is defined by the Hamilto- 
nian SHI HI 



H(t) 
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where a z = | T> <T I r 11X41 and a x = \]){[\ + ||)(T| are 
Pauli matrices. This Hamiltonian describes how at time 
t = the diabatic energies ±vt/2 of the two states cross 
with level-crossing speed v. During the crossing, the two 
diabatic states interact with a strength A, so that the 
adiabatic states (or: time-dependent eigenstates) differ 
from the diabatic states. As is usual, the adiabatic ener- 
gies are found as the time-dependent eigenvalues of the 



Hamiltonian ([T]). These are ±\/ {vt/2) 2 + A 2 , showing 
the archetypical avoided level crossing. The gap between 
these two adiabatic energies is at least 2 A, and the min- 
imum occurs at time t = 0. More intriguing is the dy- 
namics of the state 



IV(*)>=c T (t)|T> + c i (t)|i)= h§ 
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as described by the Hamiltonian {T]). Except around 
t = 0, the Hamiltonian ([T]) is dominated by the diabatic 
energies. It therefore makes sense to define an interaction 
picture by the transformation Uo(t) = exp(— \vt 2 cr z /2h), 
that is \i>{t)) = U {t)$(t)) and $(t)) = [h{t)\l) + 
£2 (£) 1 2)], where the interaction-picture probability am- 
plitudes obey 
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This system of coupled equations of motion ([3]) is one of 
the few in driven quantum mechanics that can be solved 
exactly. The dynamics at all times can be expressed in 
terms of parabolic cylinder functions [24j] . Instead of pre- 
senting this dynamics in full glory here, we will give the 
following useful summary of the dynamics 

|#=oo))=S A |V(t=-oo)), (4) 
in terms of the exact scattering matrix (or S-matrix) 



Sa = 



(5) 



Here the quantity q stands for exp(— 2irr]) and the Stokes 
phase x — tt/4 + argT(l — vq) + rj(hir] — 1), with adia- 
baticity parameter r\ = A 2 /(hv) and T the Eulcr Gamma 
function. It follows that the probability that the atom 
ends up in the initially unoccupied level ||) is given by 



P i = |c i (t = co)| 2 = l 



e ~2nA 2 /hv 



(6) 



This is the famous Landau-Zener transition probability. 
It is an exact result for all A and v. Instead of using the 
properties of parabolic cylinder functions to derive this 
result, which is the standard method, the same transi- 
tion probability for the standard Landau-Zener problem 
can be found by exact summation of an infinite series 
in time-dependent perturbation theory (38|. The latter 
method is less well known, but it is an important one, 
since it is this perturbation method that can be used to 
analyze Landau-Zener transitions in more complex sys- 
tems as well [13, Hi, Ei] , where an analysis in terms of 
special functions is not available. An example of a more 
complex system analyzed by a perturbation series is pre- 
sented below. 



III. LANDAU-ZENER TRANSITIONS OF A 
QUBIT COUPLED TO MANY OSCILLATORS 

Next we turn to our main topic, namely LZ sweeps in a 
qubit that is coupled to N cavity modes. We assume that 
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the internal interaction A of the qubit vanishes, unlike in 
the standard LZ Hamiltonian of Eq. (JTJ) , so that all bit 
flips that occur in the qubit are mediated by the cavity 
modes. The latter are described as quantum harmonic 
oscillators. The Hamiltonian becomes 



vt 
2 



N 



N 
3 = 1 



(7) 



By keeping N general here, our calculations are relevant 
both for setups with one and with two oscillators, on 
which we focus in later sections. Moreover, in Sec. lVII wc 
will show that this Hamiltonian ([7]) indeed describes the 
dynamics for LZ sweeps in circuit QED. 

We will now calculate the LZ transition probability for 
a qubit that at time t = — oo starts in its ground state 
|t). We assume that all cavity modes are initially in their 
ground states |0) as well, where |n) is a shorthand no- 
tation for the Fock state |ni,... , ri/v)- Since we chose 
A = so that the states |f) and ||) are eigenstates of the 
qubit Hamiltonian ^vta z . Then any transition between 
the two qubits states can only result from the coupling 
to the oscillator. Notably, our model |J7J is significantly 
different from the "standard" Landau-Zener problem ex- 
tended by a bath coupling via a z [H, |4l[ . 

The central quantity of interest is the probability 
Pf_^(i) that the qubit has flipped to the state ||). If 
the qubit energy is switched slowly (i.e. v — ► 0), the 
qubit will follow the adiabatic ground state which at large 
times is the state ||); then P|-_>^(oo) = 1. For large v, 
the qubit will remain in state ||), so that Pj^j(oo) = 1 
and P-|-_>.j(oo) = 0, corresponding to a nonadiabatic tran- 
sition. Generalizing the calculation presented in our 
prior work [34| from one to arbitrarily many oscilla- 
tors, we derive in the following an exact expression for 
*W°°) = EnKT,n|f/(oo,-oo)|t,0)| 2 where U(t,t') 
denotes the time-evolution operator. 

We start by a transformation to an interaction picture 
with respect to the uncoupled qubit and oscillators, 

U Q (t) = exp ( - i Vjfyjt) ex P ( - ifi vt2(J -) ( 8 ) 
i 

which yields the interaction-picture Hamiltonian 



tfW^E^l^ + ^e-^exp 



Within a perturbation expansion of the probability am- 
plitude A n = (f ,n\U(oo, — oo)|t) 0), we obtain the series 
A n = J2k==o(~ i/fi) 2ka nk- Since each H{t) flips the qubit 
exactly once, only even orders of 7 appear in A n . The 2k- 
th order contribution a n k is characterized by 2k vectors 
A^, £ = 1, ... ,2k, where A denotes a vector with exactly 
one component equal to ±1, while all other components 
vanish. The jth component {\g)j — ±1 corresponds to 



the operators b\ and bj , respectively. Then we obtain 
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2h 



^21-1) 



1=1 



(10) 



where we introduced Jl = (f2i, . . . , fijv) m order to ob- 
tain a more compact vector notation. The appearence of 
the Xi in the exponent stems from the sign in the time- 
dependent phase of the creation and annihilation oper- 
ators. The dots in the coefficient C n k(^2k, ■ ■ ■ , Ai) = 
(n\ ■ ■ ■ |0) denote the combination of 2k operators jjbj 
and -fjb^ that corresponds to the sequence A2fc, . . . , Ai. 
An important simplification of the A-summation results 
from the fact that C n k — whenever more annihilation 
than creation operators act on the oscillator ground state 
|0). Thus, we need to consider only those A-sequences 
that for all I < 2k fulfill the relation 



(ii) 



For the further evaluation, we substitute in Eq. (TT 
the times tg by the time differences Tg = tg + \ — tg, £ 
1, ... ,2k — 1, where t = t%k- Thus, we insert tg = t 
Y^t'=£ T t' so that the integral in (jTOJ) becomes 



00 oc 



2fc-l 



-00 



dt / dr 2 k-i ■ ■ ■ dri exp ft ■ Xg(t - ^ n 

" i=i e'=e 

k 2k— 1 

eXP Yl { 2T 2«-1 (* - £ T ^') - 4-1 } 



1=1 



'=2£ 



The i-integration results in the delta function 



2A- 



T 2 e-i 



E n a, 
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From the inequality (jTTJ) it follows that the second sum 
in the argument of the delta function is non-negative. 
Because the integration interval of all Tg is [0 . . . 00), any 
non-zero contribution to the integral Q12[) results from 
Ti = T3 = . . . = T2k-i = 0. Hence, the integral over 
the time differences t%, T4, . . . , Tik-2 must yield a distri- 
bution proportional to <5(ti) Sfa) ■ • • Sfak-i)- Evaluat- 
ing the integrals over all Tig separately, one finds that 
such a distribution is obtained only if Y^=i(^')j = 
for all £ — l,...,k — 1. These k — 1 relations to- 
gether with the delta function (fT3| lead to the condi- 
tions (A2f + \2t-i)j — and hence \2e-1 = —^21 for 
all I = 1, . . . , k. In combination with Eq. (fTTj) this im- 
plies that an integral is non- vanishing only if the non-zero 
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in terms of the integrated spectral density 
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(17) 
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FIG. 1: (Color online) Possible processes during LZ transi- 
tions in two-cavity circuit QED. Shown are the lowest-energy 
levels of the harmonic oscillator states \ni,n2) for TV = 2. 
The solid blue line marks a process that contributes to the 
perturbation series for Landau-Zener transition probability 
P|_,|(oo): the oscillator state jumps (repeatedly) from the 
ground state to an arbitrary one-photon state and then back 
to the ground state, in agreement with the no-go-up theo- 
rem. In contrast, the process marked by the red dashed line 
contributes to P|^f(i) only at finite times. 



The final transition probability (|16p depends only on S 
which acts as the effective coupling strength. Notice that 
quite surprisingly, the transition probability Pf^(oo) is 
independent of the oscillator frequencies Clj. Neverthe- 
less, the dynamics at intermediate times docs depend on 
them. This was shown to be the case for a qubit coupled 
to one oscillator in 34], and it will also hold for setups 
with two oscillators, as presented in Sec. fVl 



IV. MANIPULATION OF THE 
SINGLE-OSCILLATOR STATE 

We first consider LZ transitions in the standard cavity 
QED model of one qubit coupled to one oscillator. Since 
we start out in the ground state |t,0) and the Hamilto- 
nian in Eq. (0) correlates every creation or annihilation 
of a photon with a qubit flip, the resulting dynamics is 
restricted to the states |T,2n) and |J,,2n+l). As derived 
above, of the former states only |T,0) stays occupied. 
Thus, the final state exhibits a peculiar type of entangle- 
ment between the qubit and the oscillator, and can be 
written as 



+ /Pw(oc)( Cl ||l)-t- C3 |j3) + ...), 



component of \2i-1 is +1 while the same component of 
\<2t equals —1. In other words, we obtain the selection 
rule that to the occupation probability at t — 00 only 
those processes contribute in which the oscillator jumps 

(repeatedly) from the state |0) to any state with a single 

photon (i.e. to &]|0)) and back; see Fig. [TTTJ ft follows |*(oo)) = Jl - Pj^x(oo) |T0) 
that the oscillators not only start but also end in their 
ground state |0) if the final qubit state is |t). We call this 
dynamical selection rule the "no-go-up theorem" (see also 

122! )• where \c\\ 

Stating the above result in mathematical terms, of 
all possible 2fc-th order processes denoted by sequences 
of vectors -A 2 fe-i, A 2 fe-i, . . . , -A 3 , A 3 , — Ai, Ai, the pro- 
cesses that contribute to the perturbation series for 
■Pf— >T (°°) can b e characterized by simpler sequences of 
scalars jf., ■ ■ ■ , js, ji , where ji is the index of the non- 
vanishing component of both A 2 £~i and A 2 £. All prefac- 
tors C n k of contributing processes have the structure 
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1 . This implies that by measur- 
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32k-l 



2 2 jr 



(14) 



The remaining multiple integrations are performed as 
detailed in the appendix of Ref. [38[, yielding a n k = 
5n,o{Kh/v) k /k\ and is independent of the indices ji- 



Therefore the summation over ji , . . . , jk can be identi- 



fied as {J2f=i l]) k so tuat 

A n = 5 n . cxp[-TtS/hv). 
Consequently, we find the transition probability 



p w M = i-p m m = 1 



-2-irS/hv 



(15) 



(16) 



ing the qubit in state ||), a highly nonclassical oscillator 
state is produced in which only odd-photon (or: odd- 
phonon) states are occupied. Qubit and oscillator end 
up fully entangled, in the sense that after tracing out the 
oscillator states, no coherence between the qubit states 
|t) and \l) is left. The entanglement would have been 
less perfect if a nonvanishing internal interaction oc cr x 
between the qubit states |T) and |J.) had been present 
& 

While P|^|(oo) is determined by the ratio 7 2 //w, the 
coefficients C2„+i depend also on the oscillator frequency. 
In Fig. [5] we depict how for a small frequency (very small: 
equal to the coupling strength!) the average photon num- 
bers in the oscillator depend on the state of the qubit. 
In particular, the average photon number decays rapidly 
to zero in case the qubit ends "up" , in agreement with 
the no-go-up theorem derived in the previous section, 
whereas on average more than one photon resides in the 
oscillator in case the qubit has flipped to Further- 
more, it can be seen that the probability to end up in ||) 
indeed tends to the analytically exact final value. 

In contrast to the extreme situation depicted in Fig. [21 
for the recent experiments in circuit QED [l5| and for op- 
tical cavity QED one is always in the situation 7 <ti Ml, in 
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FIG. 2: (Color online) LZ dynamics of a qubit coupled to one 
oscillator, far outside the RWA regime: 7 = hQ — 0.25\/~hv. 
The red solid curve is the survival probability Pf_j(f) when 
starting in the initial state |f 0). The dotted black line is the 
exact survival probability Pj_»f (00) based on Eq. (|16p . The 
dashed purple curve depicts the average photon number in 
the oscillator if the qubit would be measured in state |j); the 
dash-dotted blue curve at the bottom shows the analogous 
average photon number in case the qubit would be measured 

IT). 



FIG. 3: (Color online) Sketch of adiabatic eigenstates during 
LZ sweep of a qubit that is coupled to one oscillator. Starting 
in the ground state |t 0) and by choosing a slow LZ sweep, a 
single photon can be created in the oscillator. Due to cavity 
decay, the 1-photon state will decay to a zero-photon state. 
Then the reverse LZ sweep creates another single photon that 
eventually decays to the initial state jj 0). This is a cycle to 
create single photons that can be repeated. 



V. ENTANGLING TWO CAVITIES BY A LZ 
SWEEP OF A QUBIT 



which case c\ w 1 to a very good approximation. Hence 
one can control via v the final state to be any superposi- 
tion of 1 10) and |J.l). In particular, in the adiabatic limit 
vh/j 2 <C 1, the final state becomes |J,1). This has the im- 
portant physical implication of the creation of exactly one 
photon in the cavity, triggered by a Landau-Zener tran- 
sition. In an experiment, the photon will subsequently 
leak out of the cavity. 

By exploiting these two processes, we propose the fol- 
lowing four-step LZ cycle for single-photon generation 
in circuit QED, as sketched in Fig. [3j The first step is 
single-photon generation in the cavity via the adiabatic 
LZ transition |f0) — > |J,1), brought about by switching 
the Josephson energy sufficiently slowly. Second, the 
photon is released from the cavity via the (controlled) 
cavity decay |J.l) — ► |J.0). In the third step, another in- 
dividual photon is generated via the reverse LZ sweep 
IJ.0) — > Ifl). Fourth and finally, a further photon decay 
completes the cycle. 

This scheme for repeated photon generation via 
Landau-Zener cycles could be implemented in circuit 
QED, where the atom-cavity coupling remains at a con- 
stant and high value and where qubits are highly tunable 
so that LZ sweeps can be made from minus to plus an 
"atomic" frequency, and back. Details about the physical 
realization in circuit QED of this proposal are deferred 
to Sec. |VlJ The strength of this scheme is its simplic- 
ity and its robustness against parameter variations, es- 
pecially variations (and fluctuations [37| ) of the oscillator 
frequency. 



Now consider the situation that the qubit is coupled 
to two cavities instead of one, with resonance frequencies 
^1,2- We will now show that the two cavity oscillators 
become entangled by a Landau-Zener sweep of the qubit, 
and moreover that the specific entangled state in which 
the oscillators end up can be engineered by varying the 
Landau-Zener sweep speed v and the frequency detuning 
8uj = (SI2 — Oi) of the oscillators. 

As before in the case of one oscillator, we will assume 
A = for the internal interaction of the qubit, so that all 
bit flips in the qubit are caused by interactions with the 
oscillators. For simplicity, we will assume that the two 
qubit-oscillator strengths are equal, 71 =72 =7- With- 
out loss of generality, we can take fli < ^2 so that the 
detuning Su> is nonnegative. The Hamiltonian becomes 

vt 

H = —a z + 7cr x (6i +b\+b 2 + b\) + hCliblh + hfl 2 blb 2 . 

(19) 

We will assume that this system starts in the ground state 
|t O1O2) before undergoing the Landau-Zener transition. 
(The subscripts 1, 2 that label the two oscillators will be 
left out below.) The general result (fT6)) then implies that 
the probability for the qubit to end up in the state ||) 
equals 

P T ^ x (oo) = 1 - Pt^tM = 1 - e ~^h 2 +f 2 )/^, (20) 

At this point it is important to appreciate that this ex- 
act result has been derived without making a rotating- 
wave approximation and by taking the full Hilbcrt space 
of the two oscillators into account. The absence of any 
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FIG. 4: (Color online) Upper panel: Adiabatic energies dur- 
ing a LZ sweep of a qubit coupled to two oscillators. Parame- 
ters: 7 = 0.25Vhv, Mli = 90Vh~v and fi 2 = WOVkv. Viewed 
on this scale of oscillator energies, the differences between ex- 
act and avoided level crossings are invisible. Lower panel: 
for the same parameters, probability Pj- -tf(t) that the sys- 
tem stays in the initial state |T00) (solid), and corresponding 
exact survival final survival probability Pf_>f(oo) of Eq. (120[) 
(dotted). 

frequency dependence in Eq. (f2"0")) is therefore quite sur- 
prising. 

In the following we are interested in the properties of 
the final qubit-two-oscillator state |V>(oo)) rather than 
merely the transition probability Pf^j(oo) of the qubit. 
In general not much can be said about this final state, but 
let us now make the realistic assumption 2 3> 7: both 
oscillator energies Mli^ are much larger than the qubit- 
oscillator couplings 71 =72 =7- Still, the frequency 
detuning Slu — (Q2 — ^1) may be larger or smaller than 
j/h. The adiabatic energies in this case are sketched in 
Fig. 21 As shown in the figure, level crossings that are 
important for the final state only occur around the times 
when the qubit energy vt is resonant with one or both 
of the oscillator energies Ml\^- There essentially only 
three qubit-oscillator states play a role in the dynamics: 
the initial no-photon state 1 | 00) and the two one-photon 
states || 10) and || 01). The most general normalized 
final state can therefore be written as 

|V(cx>)) = ^P M (co)|T00) (21) 

+ y/P^i(oo) (s 10 |110) + aoiUOl)), 

with probabilities Pf^j(oo) and Pj^(oo) given in 
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FIG. 5: (Color online) Upper panel: Adiabatic energies 
during a LZ sweep of a qubit coupled to two oscillators. 
Parameters: 7 = 0.25V~hv and O.2 = WOwHv, both as in 
Fig.Q] Hili = 80Vkv- Lower panel: Probability P T ^ T (i) that 
the system stays in the initial state | TOO) (solid), and corre- 
sponding exact survival final survival probability P^_>f (00) of 
Eq. (O (dotted). 

Eq. (|2"0"|) and with general complex coefficients s that are 
only constrained by |sio| 2 + |soi| 2 = 1 to ensure state 
normalization. Below, we will analyze the final state in 
the limits of large and vanishing detuning of the oscilla- 
tor frequencies, before addressing the intermediate case 
5u> ~ j/h. 

A. Large detuning (5u> 3> y/h) 

If the resonance energies of the cavities differ by much 
more than the qubit-oscillator coupling, then the dynam- 
ics can very well be approximated by two independent 
standard Landau-Zener transitions, see Figure \5\ The 
first transition occurs when the qubit energy vt is reso- 
nant with the energy hVli of the first oscillator. Only the 
two states || 00) and || 10) take part in this transition. 
The other transition occurs around time t — hD.%/v and 
only between the two states \\ 00) and || 01). 

This situation of two independent transitions is anal- 
ogous to optical cavity QED with two cavities and one 
atom flying successively though both of them. The 'time 
of flight' for the atom between the two cavities here cor- 
responds to the time difference St = hSuj/v between the 
two resonances. 

Just like the standard LZ transition in Sec. [HI we can 
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go to an interaction picture and summarize the transi- 
tions by S-matriccs. In the basis {|f 00), ||. 10), ||. 01)}, 
the S'-matrix for the first transition becomes 



Si= -VI 




and for the second transition 



-vT 



vY VT^7 c~ lx ~ 102 

1 

Ye ix ' + ^ 2 



(23) 

Here the quantity q' equals exp(— 2ttt]') with adiabaticity 
parameter v( = ^/ 2 /hv. The Stokes phase is now given 
by x' = + argT(l - 177') + ?/(lnr/ - 1), and fa = 
Ml^/2v. The final state after these two transitions is 

given by |^(oo)) = S2S1 — 00)). Hence, for the initial 
state \ip(— 00)) = |t 00) we find the final state after two 
transitions 



equal probabilities can not be realized. The reason is, for 
very slow LZ transitions (i.e. for q' — ► 0), all population 
already follows adiabatically the path |f 00) — + || 10) in 
the first transition, so that no population is left to take 
part in the second LZ transition and hence no population 
ends up in |J, 01). For faster transitions the situation is 
less extreme, but a population difference in the final state 
will remain. 

In summary, entangling the two cavities via a LZ 
sweep of the qubit in the case of large detuning is not 
ideal, since Bell states cannot be engineered with high 
probability because of a tradeoff between the probability 
Pf-fj^oo) = 1 — q' to create at least one photon and the 
relative probability q'/(l + q') that that one photon ends 
up in the second cavity. On the other hand, we find hat 
the state | J, 10) can be created with certainty in the adi- 
abatic limit, which means that single-photon creation in 
a single oscillator as discussed in Sec. IIVI can still be re- 
alized with LZ sweeps even if another detuned oscillator 
is present; see also |6j. 



|^(oo)) = </|T 00) 



(24) 



q'e lb 



1110) 



VII 01) 



with 8 = x' + fa- This state is a special case of Eq. (|2"Tj) 
in case of two independent LZ transitions. From the fi- 
nal three-level state (l24l) we find back the exact survival 



probability |cj(oo)|' < 



exp(— 27tr]). In the deriva- 



tion of the exact result ([20]) we took the full Hilbcrt space 
into account rather than a three-level subspace. More- 
over, we did not assume whether or not the LZ tran- 
sitions would occur independently. In other words, the 
exact probability ([20]) comes out independent of these 
two assumptions. Now that the assumptions hold, the 
exact result will still hold. However, one cannot turn 
the argument around: finding the survival probability 



|c T (oo)|< 



exp(— 2nrj) does not imply that the 



rotating-wave approximation was valid after all or that 
two LZ transitions must have occurred that were inde- 
pendent. Indeed, we will come across a counterexample 
in Sec. IVBl 

With the final state (j2"4"|) determined, one can answer 
the question whether two successive and independent LZ 
transitions are a suitable method for entangling the two 
oscillators. Conditional on measuring the qubit in state 
||), apart from an unimportant overall phase the entan- 
gled two-oscillator state becomes 



\4> osc (oc)) = 



= (110) 
+ a' V 



q'\Ql) . (25) 



In general one would like to be able to create the max- 
imally entangled symmetric Bell state (|10) + |01))/v2- 
Taking Eq. (|25[) at face value, the way to obtain equal 
probabilities for the two one-photon states would be the 
case q' = 1, but in this case the qubit would never have 
ended up in ||) in the first place, as Eq. (j2"4")l reveals. So 



B. Degenerate oscillator energies (Slo <C j/fc) 

Instead of two independent transitions, we will now 
consider the other extreme case 8u — so that during 
the LZ sweep the qubit comes into resonance with both 
oscillators at the same time. With Oi = 0,2 = the 
qubit-two-oscillator system has an extra symmetry that 
we will now exploit in the analysis of the LZ dynamics. 
Let us first go back and not yet make the rotating-wave 
approximation. The Hamiltonian (fTi?)) now becomes 



H = ycr 2 + 7^(61 + b\ +b 2 + b\) + Ml (blh + bib- 



We introduce the new operators 



b± 



V2 



(bx±b 2 ; 



(26) 



(27) 



which have standard commutation relations [6±,fr±] = 1 
and [b±,b±] = [b±,b^-] = 0. Both creation operators b± 
create one single photon with equal probability in the 
first or the second oscillator: 



bl\0 + 0-) = -j=(b\±4) |0i0 2 > 



1 

72 



(|1 1 2 )±|0 1 1 2 )). 

(28) 

So b^_ creates the symmetric and b[_ the antisymmetric 
linear combination. Instead of describing which photon 
exists in which local oscillator, we can use the fact that a 
general two-oscillator state can alternatively be written 
as 



71^=0 



t _(6j f )"+(b t _)"-|0+0-) 
V(n+)!(n_)! 



(29) 
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The above rewriting is useful because in terms of the new 
operators, the Hamiltonian (|2"6")) becomes 



hn{b\b+ + blb- 



(30) 



H = —a z + V2ja x (b + +b 

Note the factor V% in the interaction term. The Hamilto- 
nian (|30| shows that the qubit is fully decoupled from the 
antisymmetric operators. The state of the "antisymmet- 
ric photons" will therefore remain unaffected by the LZ 
sweep. Conversely, whatever the antisymmetric state the 
two oscillators are in, it will not influence the LZ dynam- 
ics. Consequently, for degenerate oscillator frequencies 
we find back the mathematical problem for a qubit cou- 
pled to one oscillator, which we already studied in [3~ij 
and in Sec. IIVI above. The difference lies in the physical 
meaning of the oscillator, either a local oscillator or a 
symmetric combination of two local oscillators. 

Now assume as before that the initial state is || O1O2) = 
|t 0+0_). Moreover, we take Mli )2 = Ml ^ 7 so that the 
rotating-wave approximation can be made. Then again 
only the three states |f O1O2), || I1O2), and || O1I2) will 
play a role in the LZ dynamics. Or, in our new repre- 
sentation, only the two states || 0+0_) and || l+0_). 
The third diabatic state || 0+l_) has a time-dependent 
energy (Ml — vt/2) and since it is annihilated by the 
RWA interaction V2j(a + b + + <r-b\), within the RWA 
the state || 0+l_) is an adiabatic eigenstate. Indeed, 
Figure [H] shows that far from the crossing time vt — 
Ml two adiabatic energies overlap with energy values 
(Ml — vt/2), which is the energy of the diabatic states 
|| I1O2) and IJ.O1I2), or equivalently, of || l+0_) and 
I J. Q_|_l_). At the crossing, one such a diabatic line re- 
mains and it corresponds to the state || 0+l_). The 
other two states |f 0+0-) and || 1+0-) form the avoided 
crossing in Fig. [6] 

The LZ transition probabilities can again be calculated 
with the help of an 5-matrix. The (2 x 2)-dimensional 
S-matrix in the {|f 0+0_), || l+0_)} basis gets the form 



q" e~ lx 



q" e 1 * 



(31) 



where q" = exp(— 2nr]") and the Stokes phase is given by 
X" = 7r/4 + argT(l-ir/')+?7"(ln77"-l) + Ml 2 /2w, with 
adiabaticity parameter 77" — (V2"/) 2 /hv. It immediately 
follows that |c|(oo)| 2 = q" = exp(— ATrj 2 /hv), again in 
agreement with the exact result of Eq. ([20]) . Thereby 
the same generally valid exact result indeed shows up in 
the two limiting cases that can be treated with simple 
S-matrices: the case (Mli j2 7, 8u) 7) in Sec. IV Al 
and here the case (Mli = Ml 2 S> 7). 

The final three-level state in these two limiting cases is 
different, though. Instead of the state (f2"4")) . for the qubit 
coupled to degenerate oscillators starting in || O1O2) we 
now find the final state 



|V(oc)) = v^lt 0+0_) - y/T^qtjx ||l+0_) 
= v^lt 0x0 2 ) - y/l=?e«' 
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FIG. 6: (Color online) Upper panel: Adiabatic energies 
during a LZ sweep of a qubit coupled to two oscillators 
with degenerate energies. Parameters: 7 = 0.25V^v and 
MI2 — lOOVm;, as before; this time Mli = hQ.2- Lower panel: 
Probability Pf _>-)■(£) that the system stays in the initial state 
I L00) (solid), and corresponding exact survival final survival 
probability P T _ T (oo) of Eq. ((20j) (dotted). 



Clearly, if the qubit is finally measured and found in state 
||), then the two oscillators end up in the symmetric Bell 
state 

|^ osc (oo)) = -^(|li0 2 ) + |0il 2 )). (33) 



Notice that unlike in the final oscillator state ([23]) af- 
ter two independent LZ transitions, no relative phase be- 
tween the two oscillator states | I1O2) and |0i I2) is built 
up here. For degenerate oscillators this is as one would 
expect. 

It follows from Eq. (f3"2"]) that in practice one has two 
options to produce this Bell state (f3"3"]) : The first option 
is to employ fast LZ transitions which produce either the 
sought oscillator Bell state or the oscillator ground state. 
A measurement of the qubit state is then required, as in 
a proposal {Toj to create Bell states in optical two-cavity 
QED. If |t) is measured, then a fast LZ sweep back to 
the initial ground state is needed and the process can be 
repeated until the qubit is measured in state ||), whereby 
the Bell state is produced. The second option is simpler: 
one chooses slow LZ transitions and by adiabatic follow- 
ing produces the Bell state (j33|) with almost certainty. 
Me reason why the first option may be preferred after 
til in practice is that the usual fight against decoherence 
iay require faster than adiabatic operations. 
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FIG. 7: (Color online) Upper panel: Adiabatic energies 
during a LZ sweep of a qubit coupled to two oscillators 
with large energies, and with detunings of the order of the 
qubit-oscillator coupling 7. Parameters: 7 = 0.25^/hv and 
— 100\/fiu, as before; M7i — 96\/fru. Lower panel: Prob- 
ability Pf_>f (t) that the system stays in the initial state | f 00) 
(solid), and corresponding exact survival final survival prob- 
ability P T _ T (oo) of Eq. (HI) (dotted). 



An important difference with the final oscillator 
state (I25p after two independent LZ transitions in 
Sec. IV Al is that in degenerate oscillators always the sym- 
metric Bell state will be produced by measuring the qubit 
state in ||), independent of the adiabaticity parameter 
77". So if the goal is to produce the symmetric Bell state, 
then one should try and build systems where a qubit 
couples with equal strength 7 to two oscillators with fre- 
quency detuning Slu <C 7. If on the other hand one would 
like to be able to vary the final two-oscillator state, then 
detuned oscillators are to be preferred that give two in- 
dependent LZ transitions, producing the state ([^5]) . 



C. Intermediate detuning (Slu ~ 7/ft) 

In case of intermediate detuning, the dynamics cannot 
be described as one or as two independent successive LZ 
transitions (see Fig. [7]). Of course we know from the exact 
result that the survival Pj_^f (00) probability is given by 
Eq. (f20|) . as before. Apart from this probability, we are 
again interested in the final state. If we assume again that 
^1,2 3> 7 and make the rotating-wave approximation, 
in the local-oscillator basis we end up with the (3 x 3) 



Hamiltonian 

(vt/2 7 7 \ 

H RWA (t) = 7 m- vt/2 . (34) 

\ 7 nn 2 - vt/2 J 

This is a special case of the Demkov-Osherov model [42j , 
in which one level crosses N parallel levels. Interestingly, 
the transition probabilities (fTB]) and (l2"0|) are also exact 
within the Demkov-Osherov model [42|, |43j. Even the 
scattering matrix S for the (iV+1) levels is known exactly 
[44I [4{| . The interesting result is that the final state is 
still given by Eq. , as if the two level crossings had 
been independent. Only the phase 9 in Eq. is to 
be replaced by a more complicated expression (45|, but 
for the two-oscillator state |V> sc(oo)) this is an irrelevant 
overall phase. A reason to avoid intermediate detunings 
in experiments is that convergence of the relative phase 
to the final value (4>2 — 4>i) in Eq. ([24")) is reached more 
slowly than for large detuning [4J|, whereas (02 — 4>i) 
simply vanishes for zero detuning. 

In these Sees. IV AllV CI we considered three regimes for 
the detuning, but we have not yet estimated for which 
parameters these regimes occur. For example, in practice 
the two oscillator frequencies will never be exactly degen- 
erate, so the question arises how to define the regime in 
which they are effectively degenerate. This can be es- 
timated by requiring that the typical time of a single 
LZ transition tlz = 2j/v [H, [4(| is much smaller than 
the sweeping time HScj/v from energy Ml\ to Ml2- The 
sweep velocity drops out and we find the requirement 
5u> <C 2^/h. Following a similar reasoning, the regime 
of two independent LZ transitions is characterized by 
Slu > 27/fi. 

VI. EXPERIMENTAL REALISATION 

Above we stated that the Hamiltonian ([7]) can be re- 
alized experimentally, and here we briefly outline a real- 
ization in circuit QED. Several realizations are promising 
0, EE EH E3, EH ' an d nere we consider in more detail 
the setup of the experiments at Yale [HI, [lj| • 

We first consider the setup with one qubit and one 
oscillator, which has already been realized and for which 
we proposed in Ref. and in Sec. HVI to generate single 
photons with LZ sweeps. In tunneling representation, a 
circuit QED setup with one circuit oscillator is described 
by the Hamiltonian (l5l [l6j 

H(t) = -^a x -^p-a z +hntfb +1 (tf+b) [a x + l- 2N g ] 

(35) 

for the qubit, the circuit oscillator, and their mutual 
coupling. The circuit cavity can indeed be modelled as 
a harmonic oscillator. The electrostatic energy E c \ = 
4E C [1 — 2N g ] is determined by the charging energy E c 
and the tunable gate charge N g . The tunable flux <I>(i) 
penetrating the superconducting loop will be used to 
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drive the qubit. The flux controls the Josephson energy 
Ej(t) = Sj )TOax cos[7r$(t)/$o]) where <& is the flux quan- 
tum. The two-level approximation underlying the Hamil- 
tonian (|35[) is valid in the "charge regime" E c ^> Ej. 
In order to minimize dccoherence, one typically oper- 
ates the qubit at the optimal working point N g = 4, 
so that E c i — [47|- Here we also restrict ourselves to 
this optimal working point. The LZ dynamics can then 
be realised by switching the flux <I>(i) in such a way that 
Ej(t) = —vt, with v > 0. This way, our central Hamil- 
tonian ([7]) is realized for one qubit coupled to one oscil- 
lator (TV = 1). Moreover, the very low temperatures in 
circuit QED experiments [l5| justify the assumption for 
our calculations that both the qubit and the oscillator 
are initially in their ground states, i.e. \^(— oo)) = |f, 0), 
where a z \~\) = |T>- 

For an ideal LZ sweep, the Josephson energy should 
be swept from Ej — — oo to Ej — oo. In practice, 
the two-level approximation is valid in a finite energy 
range. Moreover, in reality Ej is bounded by an Ej iiaax 
which is determined by the critical current. The condi- 
tion -Ej imax > Ml is required so that the qubit comes into 
resonance with the oscillator sometime during the sweep. 
The duration of this linear sweep has to be long enough, 
so that transition probabilities have converged and the 
finite time interval can be extended to t — — oo . . . oo 
in calculations describing the dynamics. In circuit QED 
this situation occurs, since qubit energies can be swept 
around the oscillator resonance Ml over intervals much 
larger than the interaction 7. 

As another practical condition, inverting the flux 
through the superconducting loop requires a finite time 
2T min , so that v cannot exceed v max = E^ max /2T min . For 
the setup of Refs. [H, EH, the sign of the initial Joseph- 
son energy -Ej !max = 2nhx 10 10 Hz can be inverted within 
T = 1/xs so that u max = 2ith x 10 16 s -2 . 

The cavity frequency fl and the qubit-oscillator cou- 
pling 7 are determined by the design of the setup. A 
typical cavity frequency is f2 = 2tt x 10 9 Hz. For the 
qubit-oscillator coupling strength we assume ~//2irh = 
3 x 10 6 Hz. Since 7 -C Ml, the generation of more than 
one photon is negligible. In this limit, by choosing a 
proper value of v, one can obtain any desired superposi- 
tion of the states |f0) and ||1), as explored in more detail 
in H|. 

In Sec. IV Bl we proposed to create Bell states by cou- 
pling the qubit with identical couplings 7 to two degener- 
ate noninteracting circuit oscillators. Such a setup with 
one qubit and two oscillators currently does not exist, 
but a detailed proposal how one could fabricate one is 
given in Ref. [21(. Our wishlist of (almost) identical cou- 
plings 71,2 and identical resonance frequencies flip may 
not be realized easily. The frequencies are typically are 
at least 10 4 times larger than the couplings. Ideally one 
of the two oscillator frequencies is tunable, so that it can 
be brought into resonance with the other one. 



VII. SUMMARY AND CONCLUSIONS 

We have calculated the exact Landau-Zener transition 
probability of a qubit that is coupled to one or to sev- 
eral quantum harmonic oscillators, given that the system 
starts in the ground state. The resulting transition prob- 
ability does not depend on the oscillator frequencies. 

The LZ dynamics of a qubit coupled to one oscillator 
can therefore be manipulated via the sweep speed of the 
transition. In state-of-the art circuit QED, both adia- 
batic and nonadiabatic transitions could be realized ex- 
perimentally by varying the magnetic flux. For different 
speeds of the transition, the qubit and the oscillator end 
up in different entangled states, so that LZ transitions 
can be part of the toolbox to prepare qubit-oscillator en- 
tangled states. 

Moreover, we have shown how LZ transitions in a qubit 
can be employed to entangle two oscillators. Especially 
for oscillators with equal energies, it was found that LZ 
transitions are a robust way to create in the oscillators 
the maximally entangled state (|01) + \10))/y/2 that is 
known as the symmetric Bell state. The generation of 
this particular state requires the qubit to undergo a LZ 
transition induced by both oscillators at the same time. 
In circuit QED such a situation may be engineered, for 
example by coupling two orthogonal transmission-line 
resonators to the same superconducting qubit [2lj . 

In standard cavity QED, maximally entangled cavity 
states have been realized experimentally, but only for two 
degenerate optical modes of the same optical cavity 
Our proposal is different in that we consider two spatially 
separated circuit oscillators, which is somewhat analo- 
gous to the optical protocol proposed in Refs. [13, Ell ]. 
The difference between the present work and another re- 
cent proposal to create Bell states in superconducting 
circuits [481 ] is that we consider entanglement creation in 
non-interacting cavities, whereas in [48| it is shown how 
Bell states could be created in capacitively interacting 
qubits. 

The decisive advantage of our proposal is that qubit- 
oscillator interaction strengths are static. The archetyp- 
ical way to entangle two optical cavities is by passing an 
atom successively through both cavities. This leads to 
one (non-)adiabatic transition followed by another. Ex- 
perimentally, there will be a spread in the velocity of the 
atom when repeating the experiment. The analogue of 
this situation in circuit QED that we have studied here 
does not suffer from this disadvantage: a single qubit 
can be swept through two oscillator resonances either si- 
multaneously or successively, in both cases with a well- 
controlled constant sweep speed. 

In principle, our proposal to couple two circuit os- 
cillators to a qubit to create Bell states in them can 
be generalized to three or more circuit oscillators. For 
three frequency-degenerate oscillators with equal cou- 
pling strength 7 to the qubit, instead of the symmetric 
Bell state after the LZ sweep of the qubit, one would 
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obtain the so-called W state [49|, 15CJ 

\W) = (|001) + |010) + |100» /\/3, (36) 

while for N oscillators the iV-qubit W state would be 
produced [13, IKH ■ Needless to say that it is techno- 
logically very challenging to fabricate such systems, but 
the creation of these maximally entangled states simply 
by a single LZ sweep in a qubit would be fascinating. 
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